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The self consistent effective medium approximation (SEMA) is used to study three-dimensional 
random conducting composites under the influence of a strong magnetic field B, in the case where 
all constituents exhibit isotropic response. Asymptotic analysis is used to obtain almost closed 
form results for the strong field magnetoresistance and Hall resistance in various types of two- 
■ and three-constituent isotropic mixtures for the entire range of compositions. Numerical solutions 

of the SEMA equations are also obtained, in some cases, and compared with those results. In 
two-constituent free-electron- metal/perfect-insulator mixtures, the magnetoresistance is asymptot- 
U ' ically proportional to |B| at all concentrations above the percolation threshold. In three-constituent 

metal/insulator/superconductor mixtures a line of critical points is found, where the strong field 
magnetoresistance switches abruptly from saturating to non-saturating dependence on |B|, at a 
certain value of the insulator-to-superconductor concentration ratio. This transition appears to be 
related to the phenomenon of anisotropic percolation. 

To appear in Phys. Rev. B 
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Over the years, the Bruggeman self consistent effective medium approximation (SEMA) for the bulk effective 
electrical conductivity of a composite mediumllm has spawned a number of extensions and generalizations. Those 
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include a SEMA for elastic stiffness of a composite mediumPtS a SEMA for_a-.weakly nonlinear conducting composite,^ 
and a SEMA for a strongly nonlinear (power law) conducting composite.Qu They also include a SEMA for.-a-,linear 
conducting composite medium where the constituents are characterized by nonscalar conductivity tensorsaxi. £Luch 
approximations are needed to study the nonscalar conductivity due to an externally applied magnetic field B.ETllj 
' In composite conductors, the microstructure can induce a dependence of the Ohmic resistivity on magnetic field 
(i.e., magnetoresistance) even when none of the constituents exhibits any such dependence by itself. This arises 
because the local current density J(r) acquires spatial fluctuations in both magnitude and direction, as a result of 
the heterogeneity or microstructure. Those are reflected in similar fluctuations of the local Hall (electric) field E#(r), 
which usually has a nonzero component along the direction of the volume averaged current density (J). The volume 
average of that component of E# vanishes to leading order in B, but in higher orders that average is usually nonzero. 
As a result, the bulk effective Ohmic resistivity of the composite usually depends upon B. The SEMA developed 
g ' in Refs. |],[| was used in numerical studies o^p^gne^otransport in various types of composites, where this induced 



magnetoresistance played an important roleJaH'Ejfcj It was also used in an unpublished discussion p£ asymptotic 
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■ strong field behavior of magnetoresistance in a free-electron-conductor host with open orbit inclusions.! 

More recently, it was shown that this type of SEMA, when applied to composites with a columnar microstructure, 
often violates some exact relations which exist between the different components of the bulk effective resistivity 
tensor. A modified SEMA has been developed which incorporates those exact relations [the "columnar unambiguous 
self-consistent effective medium approximation (CUSEMA)"]; it has been used for a detailed asymptotic analysis 
of metal/insulator (M/I) and metal/superconductor (M/S) random columnar mixtures in the strong field limit.Ej 
Very recently, the CUSEMA was applied to a three-constituent metal/insulator/superconductor (M/I / S) random 
columnar composite mixture; a line of critical points was found to appear whenever the / and the 5* constituents are 
present in equal amounts, i.e., when their volume fractions pi and ps are equal. At this point, the asymptotic strong 
field behavioc-of the magnetoresistance switches abruptly from saturated for ps > pi to nonsaturating oc |B| 2 behavior 
for ps < p/.c-j This critical point exhibits scaling behavior, along with critical exponents and scaling functions, all of 
which were determined, approximately, using CUSEMA. 

In this article we present a theoretical study of magnetoresistance in three-dimensional disordered composite media. 
Some of the systems we study are isotropic, two-constituent mixtures of isotropic constituents. In this category, we 
consider mixtures of normal conductors (M1/M2), M/I mixtures, and M/S mixtures. No "intrinsic anisotropy" is 
allowed for any constituent; thus, open orbit constituents arc not considered. The normal conducting constituents 
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can be simple free-electron conductors, in which case only the Hall resistivity depends upon B, or they can be more 
complicated conductors or semiconductors where the Ohmic resistivity also depends upon B. But the longitudinal 
Ohmic resistivity pn (along B) and the transverse Ohmic resistivity p± (perpendicular to B), as well as the Hall 
resistivity pHaii, are assumed to be independent of the direction of B. We also consider an M/I / S mixture, again with 
isotropic constituents, an isotropic microstructure, and Py > P±\ Puali °^ ^ ne ^ cons tituent which are independent 
of the direction of B. 

The SEMA equations for these composites do not admit exact closed form solutions. However, in the strong field 
limit, asymptotic analysis can be applied to those equations. The strong field limit means that \H\ 3> 1, where H is 
the Hall-to-transverse-Ohmic resistivity ratio in a metallic or normal conducting constituent 

(M) 

H=J$- = p\B\, (1.1) 

Pi 

and where p is the Hall mobility; note that p and p^fn, and therefore also H, can be either positive or negative, 
depending on the sign of the majority charge carriers. This asymptotic analysis often leads to results in simple closed 
form, or at least "almost closed form." In several cases we have also solved the SEMA equations numerically in order 
to compare with the asymptotic analysis. 

The rest of this article is organized as follows. In Section [n] we briefly review the SEMA, using a formulation that 
reproduces the usual SEMA equations for non-scalar conducting constituents in terms of depolarization coefficients 



of the different types of inclusions. In Sections III and OV asymptotic analysis is used to find almost closed form 



solutions to those equations for a number of special cases. Those asymptotic solutions are compared, in some cases, 



with numerical solutions of the same equations. In Section [II, we treat the special case of an M1/M2 binary mixture 



All concentrations are considered, starting from the limit of extreme dilution and all the way to the percolation 



threshold. In Section IV, M/I / S three-constituent mixtures are treated using asymptotic analysis. A line of critical 
points is found where the strong field response changes abruptly in the manner of a second-order-phase-transition, as 
the relative amounts of / and S are varied. Critical exponents, a scaling variable, and scaling functions are found which 
characterize the critical behavior. We also discuss M/S and M/I mixtures as special cases of these three-constituent 
composites. A linear dependence of the magnetoresistance upon |B| is found in free-electron-metal/perfect-insulator 
mixtures for all concentrations where the mixture conducts macroscopically. 

Section [v] provides a discussion of the results. We formulate a physical picture of the microscopic processes (i.e., 
local current flow patterns) which lead to some of the results found in previous sections for the macroscopic response. 
We also discuss a possible relation between the line of critical points found in the macroscopic magnetotransport of 
three-constituent M/I/S mixtures and the phenomenon of anisotropic percolation. 



II. REVIEW OF SEMA FOR THE CURRENT PROBLEMS 



The self consistent effective medium approximation (SEMA) for constituents with arbitrary conductivity tensors, 
which appear in the system as ellipsoidal grains, was developed many years ago by one_pf the present authorsEI. 
This development-, followed more specialized versions, such as all scalar conductivities pa and later, strong field 
magnetotransportfi and weak field Hall conductivity.!!!] Here we describe a slightly different formulation of the general 
theory of Ref. ^. We use this formulation to find asymptotic physical solutions of the SEMA equations when the 
magnetic field is very strong. This formulation was described in detail in Ref. [l6|. A similar formulation also appeared 
in Ref. [if 

In practice, the SEMA requires that one calculate the electric field Ei induced inside a single inclusion, with 
conductivity tensor <7i nc , embedded in an otherwise uniform host with conductivity tensor <Thost, when an external 
uniform electric field Eo is applied at large distances. Whenever the inclusion is an ellipsoid, Ei is uniform, whatever 
the values of a- iric , (Thost- Obviously, Ei will be a homogeneous linear function of Eo, which can be written with the 
help of a matrix 7 in c(<Tinc, Sw) 

E^TJnc-Eo. (2.1) 

If the coordinate axes are taken to lie along the principal axes of the symmetric part of the tensor Ohost, then we can 
write 

1 ^ =*~r-E , n ^ (2.2) 

cry a ( _(host) (host) \ 

P \Uaa Q8 ) 

5& = (Thost — (Tine- (2.3) 
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The factors n a p are elements of the depolarization tensor h, not of the actual physical inclusion, but of its image after 
the coordinate axes have been rescaled by the following transformation 



(host) 
0~ xx 



y 



(host) 

a yy 



(host) 



(2.4) 



Under this transformation, the ellipsoidal inclusion is usually transformed into another ellipsoidal shape, with major 
axes that have different lengths and different orientations compared to those of the physical inclusion. The SEMA 
is then obtained by setting <7h os t = &ei an d demanding that the small but exactly calculable change in (J), caused 
by one isolated ellipsoidal inclusion, vanish when averaged over the different kinds of inclusions. This leads to the 
following self consistency equations for the elements of a e 
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• 7inc(o"inci & e)) 



(2.5) 



In general, these are nonlinear equations, in which there appear non- elementary transcendental functions of the 
elements of a e . Thus, closed form solutions are usually out of the question. However, we will show below that 
asymptotic solutions, for a very strong magnetic field B, can sometimes be obtained in almost closed form. For 
smaller values of B, the SEMA equations can be solved numerically, where needed; these solutions are discussed 
below. 

We assume the following forms for the resistivity matrices of the host and inclusion: 



Phost 





(2.6) 



These forms mean that both the host and the inclusion have isotropic electrical response, and that the only physi- 
cally selected direction is that of the magnetic field B, which is always taken to lie along z. This implies that the 
microstructure, as well as the electrical response of all the constituents, are isotropic. The a coefficient in phost is 
actually redundant. We will sometimes exploit that redundancy by choosing po of the bulk effective resistivity tensor 
of the composite medium to be the same as one of the constituent values of p\. The conductivity tensors cxhost, 0mc 
are easily found by inverting phost, Pine- 

Because we are considering isotropic microstructures, we will assume that all the constituents appear as spherical 
inclusions in the fictitious uniform effective medium host. In that case we find 



7inc(Omc, Chost ) 
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where 
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and where n x , n y — n x , an d n, are the depolarization factors of the spheroidal shape into which the spherical inclusion 
was transformed by Eqs. j2.4| ). Those depolarization factors are elementary transcendental functions of a, (3, A.EJ 
Thus, the elements of 7i nc are also non-algebraic functions of those parameters. Con sequ ently, the equations obtained 
for those parameters, by applying the self consistency requirements implied by Eq. (2.5), will usually be complicated, 
coupled, non-algebraic equations which lack closed form solutions in terms of elementary functions. This makes the 
qualitative study of their physical solutions highly nontrivial. 



III. APPLICATION OF SEMA TO MAGNETORESISTANCE OF AN Mi /Mi MIXTURE 

As an illustration, we now work out the effective resistivity tensor for a binary composite of two normal metals, 
using the SEMA. Two other special cases — a composite of normal metal and insulator, and one of normal metal and 
perfect conductor — will be discussed in the next section as special cases of three-constituent mixtures. 

In the present case, the bulk effective resistivity tensors of the two constituents are assumed to have the similar 
forms 
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P2 = P2 
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while the bulk effective resistivity tensor is assumed to have the form 



(3.1) 




(3.2) 



Note that p e and p\ have been expressed in terms of the same resistivity factor p\. The forms (3.1) and (3.2) imply 
that both constituents are isotropic conductors, as is the composite, and that the magnetic field B is applied along z. 
They include the case where both p\ and p 2 represent simple free-electron or free-hole conductors. They also include 
more general types of con ductors, where the Ohmic resistivities and Hall mobility depend upon |B|. 

Application of Eq. ( |2.5| ) leads to three coupled equations for the unknown parameters a, 0, A, which arise from the 
zz, xy, and xx components of that tensorial equation [the xx and yy comp one nts lead to the same equation, as do the 
xy and yx components, while the xz, zx, yz, and zy components of Eq. (2.5) vanish identically; p\ and p 2 = 1 — Pi 
denote the volume fractions of the two constituents] : 
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Not surprisingly, these equations do not admit any closed form solutions. Since we are interested in the behavior of 
a, (3, A at strong fields, we will apply asymptotic analysis for the limiting case where Hi = H and H 2 = luH are both 
much greater than 1, i.e., ui remains bounded but \H \ 3> 1. All other physical parameters of the two constituents 
are also assumed to remain bounded. From Eq. (3.3) it follows that A must also remain bounded. We have tried 
various ansatzes for the asymptotic behavior of a, f3, and proceeded to examine whether they lead to a mathematically 
consistent and physically admissible solution. The only ansatz that could satisfy those requirements was 



P^poH, a = ao\H\ 2 ^, A = A , 



(3.6) 



where cto, Po, Xq are positive dimensionless coefficients. It follows from this ansatz that the transformation of Eqs. 
( |2.4p changes the spherical inclusions into flat, pancake shaped or oblate spheroids, with the short major axis iStlpng 
z, and with eccentricity e that diverges as \H\ — > 00. The eccentricity e and depolarization factors are given bylia 
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(3.7) 
(3.8) 
(3.9) 
(3.10) 
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This leads to D\ = 1 and D2 = 1, and when these results are used in Eqs. (3.4) and ( |3.5| ), one arrives at the following 
explicit asymptotic expressions for ao, (3q, Xq: 



/3 = 



1 



A =Pl^H P2^2, 

Pi 



P_2_ 

0\ 
2/3 



(3.11) 
(3.12) 
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(3.13) 



Note that ao = when the Hall resistivities of the two constituents are equal U\p\ = Hipi- Moreover, when 
Pi I ' p\ 1, then /3o tends to the finite value 1/pi, while Ao and ao both diverge. As we shall see in Section [V C 2 
below, this behavior presages the behavior exhibited by those parameters when the No. 2 constituent is a perfect 
insulator, i.e., when P2 = 00. 

A similar asymptotic dependence of the transverse relative bulk effective magnetoresistance on magnetic field was 
found previously for a free-electron- metal host with open orbit inclusions .113 A similar asymptotic dependence was also 
previously found by Dreizin and Dykhne, who presented a qualitative microscopic discussion of magnetotransport in 
composite media. t3 In our language, those results would translate into a oc \H\ 2 ^ 3 , as we have also found. While the 
asymptotic behavior found in those previous studies is similar to what we find here, the present discussion shows that 
(a) one does not need to have open orbits in any constituent in order to observe this kind of behavior, (b) the behavior 
we find is a straightforward consequence of SEMA. Since SEMA becomes exact in the dilute limit, when either pi 
or P2 is very small, we believe that the asymptotic result a oc |iJ| 2 / 3 is exact. It appears that the only requirement 
for obtaining this type of asymptotic response is that the two constituents have comparable Ohmic resistivities and 
different Hall resistivities. In fact, our numerical solutions of these equations also suggest that only the difference in 
the Hall resistivities is crucial: The |_ff| 2 / 3 power law is obtained even if the two constituents have the same Ohmic 
resistivities. 

In order to confirm that our ansatz for the asymptotic behavior is indeed correct, we have solved the SEMA 
equations numerically for a simple example. We assume that pi = p%, v\ = v% = 1, and H = Hi = 2H2, and 
we consider pi = P2 = 0.5. In Fig. |l](a) we plot the relative bulk effective transverse magnetoresistivity a — 1 vs. 
|i?| 2 / 3 . It is evident that this quantity rapidly approaches a linear dependence on |i?| 2 / 3 as \H\ increases, and that the 
asymptotic dependence found analytically above is accurately achieved for \H\ > 2 for this choice of parameters. In 

Fig. |l|(b), we plot the quotient 0/H = ^Haii/PHaii vs ' \H\ 2 ^ 3 - The second form of this quotient shows that it is equal to 
the ratio of bulk effective Hall resistivity to the Hall resistivity of the No. 1 constituent, suggesting the name "relative 
Hall resistivity" . Clearly, (3 also rapidly approaches its asymptotic dependence upon H , which is linear rather than 
oc |£f| 2 / 3 i.e., the ratio (3/H becomes field- independent. 



IV. THREE-CONSTITUENT M/I/S MIXTURES 
A. General considerations 



Next, we turn to a discussion of a three-constituent composite containing a volume fraction pu of normal metal, 
Ps of perfect conductor, and pi of insulator. As part of this discussion, we will consider the special cases of two- 
constituent N/S and N/I mixtures. We use the subscripts M, I, and 5* to denote normal metal, insulator, and perfect 
conductor. 

We assume the following forms and inequalities for the bulk effective resistivity tensor p e , and for the three con- 
stituent resistivity tensors /5m, Pi, Ps'- 
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Pm = Pm 




Ps 



= p s i, pi = pil, ps < Pm < \H\p M < Pi, 



(4.1) 
(4.2) 
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where 1 is the unit tensor. We use the results of Eqs. (2.7) and (2.8) in Eq. (2.5) for this system, and then take the 
limits ps — > and pi — > oo to get the followi ng e quations for a, P, A from the zz, xy, and xx components of that 
equation [again, the other components of Eq. (2J3) lead either to redundancies or to identities; pm, Pi, Ps denote the 
constituent volume fractions, which satisfy pu + Pi + Ps = 1] 
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Combining Eqs. (^J) and ( [4.5|) we can simplify Eq. (|4.5| ) slightly to get 

i 2 + f3 2 ( H\ n 2(3H 
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Despite the seemingly simple form of Eqs. ( [4.3j ) and (4.4), the right hand sides of those equations do not constitute 



explicit expression s for A or for P/(a 2 + /3 2 ), because those right hand sides depend upon a, P, A through n z and 



Nevertheless, Eq. (4.3) does show that n z must satisfy the following inequalities in order for A to be positive: 



Ps < n z < p s + pm — 1 - Pi, 



(4.9) 



where the lower bound must satisfy ps < p c = 1/3 and the upper bound must satisfy ps + Pm > Pc = 1/3 (as before, 
p c denotes the percolation threshold, equal to 1/3 in SEMA) in order for the entire composite to have a finite, nonzero 
bulk effective conductivity. 

Eqs. ( |4.3| )-(4.8) do not admit closed form solutions. In order to obtain asymptotic solutions when \H\ ^> 1, we tried 
a range of possible ansatzes for the asymptotic forms of a, 0, A, and proceeded to examine them for mathematical 
consistency and physical admissibility. Only three of those ansatzes stand up to both requirements; each is valid for a 
different range of the constituent volume fractions pu, Pi, Ps, where, of course, pu +Pi +Ps = 1- The three ansatzes 
and the resulting solutions are described in the following subsections. 



B. The saturating regime 

1. General case 



The first ansatz that leads to admissible results is 



a = a , A = A , 



Po 
H' 



(4.10) 



This leads to the following results: 
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D M = (1 - n x f, Di S (1 - n x )* 4^511 
Consequently, Eq. (4.4) yields the following relation between /3q and ao 

a 2 PM 

Pa = a Q ■ . 

Pi + Pm 



(4.11) 



(4.12) 



Using this to eliminate @ from Eq. (4.8), we get 

1 - n 
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^> n z = 1- 2n x = 1 - 2^5 = p M +PI -ps > 



(4.13) 



A consequence of the last inequality for n z is that the transformed shape of the spherical inclusions is an oblate 
spheroid. Since we must also have [see Eq. (4.9)] 



1 - 2p s = n z <ps + PM = I- Pi, 



(4.14) 



we find that we need to require 2p$ > pi in order for this asymptotic solution to be valid. Eq. (4.3) now becomes 

2p s (l - 3ps) 



Ao = v 



{1 - 2p s ){2ps ~ pi)' 



(4.15) 



In order to determine ao we first need to solve the following transcendental equation for the eccentricity e of the 
oblate spheroid: 



1 + e 2 1 

5 — (e — arctane) = 1 — 2ps > 

e J 3 



(4.16) 



and then use the relation between e and a, /3, X 
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a = A (l + eg). 

These asymptotic results are valid for the range of constituent compositions defined by 

Pi 1 

T <P5< 3- 

When ps pi/2 from above, ao, Ao, and (3q all diverge, but at different rates: 
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2ps - Pi 
1 

(2p s -pi) 2 ' 
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(4.17) 
(4.18) 

(4.19) 
(4.20) 



It is worth noting that a and A are both proportional to v. Along with the fact that they are asymptotically 
independent of H, this indicates that the local electric field in the M constituent is directed mostly along B. Similar 
proportionality is also found in the two-constituent M / S case, discussed below. It is also worth noting that when 
Ps — > 1/3 from below, then also n z — * 1/3, and consequently eo — * 0. Both ao and Ao tend to as 1 — 3ps, but 
oo/Ao — * 1- This indicates that when the S constituent approaches its percolation threshold p c = 1/3, the same 
current flow paths in the M constituent are responsible for the leading contribution to the macroscopic response 
whatever the direction of the average current density (J) . 
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2. Two-constituent M/S mixture 



An important special case of a three-constituent composite with saturating behavior is pi — 0, cor resp o ndin g to 



a two-constituent M/S mixture. One can also get t he re sults for this case by setting pj = in Eqs. (f4.3|), (4.4) [or 



alternatively by taking the limit pxj pi — > oo in Eqs. (3.3), ( p.4| )1 to obtain the following equations: 



and noting that Eq. (4.5) [or Eq. (3.5)] is satisfied as an identity to leading order in 1/H 

The asymptotic large \H\ behavior can then be obtained by using the ansatz of Eqs. ( |4.10| ). The second of Eqs. 



(4.21) immediately leads to a simple relation between /3o and ao, namely 

A) - ag. (4.22) 

The value of the ratio ao/Ao is obtained by first solving the transcendental equation for the asymptotic eccentricity 
eo, namely 

1 + e 2 

n z = p M ~Ps = 5-^(eo - arctaneo), (4.23) 



and then using 

Ao 



O "- (4-24) 



Finally the first of Eqs. (4.21) leads to 



A -,(l-— ^ =i J-^L. (4.25) 

V PM - PS/ PM ~ PS 

The fact that ao, /3o, Ao all depend upon i.e., ao oc Ao oc (3q oc v 2 , indicates that the leading contribution to 
the macroscopic response is due to local currents that flow parallel to B in the M constituent. 

When ps approaches the SEMA percolation threshold value 1/3 from below, then n z — > 1/3; hence eo — > and 
ao/Ao — ► 1 from above. Both ao and Ao tend to as 1 — 3ps, while (3q — > as (1 — 3ps) 2 - The fact that Ao/ao — * 1 
again indicates that the same current flow paths in the M constituent are responsible for the leading contribution to 
the macroscopic response whatever the direction of the average current density (J)|— ■ 

The SEMA equations for an M/S mixture were treated numerically in the past.O Those calculations are in agree- 
ment with the asymptotic pm dependence of ao, Ao, and /3o obtained for such mixtures in this subsection. 

C. The non-saturating regime 

1. General case 

Another ansatz which leads to an admissible solution is 

a^a \H\, X^Xo\H\, (3 (3 H. (4.26) 



Because A is now very large, Eq. (4.3) leads to the following results 



n t *p s +PM-o(±^ = i- P j-o(j f A , (4.27) 

"■^-^-T + Kilf)- (4 ' 28) 

Since pj < 2/3, we must have n z > 1/3 if \H\ is large enough, and the transformed spherical inclusion is then again 
an oblate spheroid. 

Instead of using the unknowns ao, /3q, it is convenient to introduce the variables x and y, defined by 
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Eq. (4.4) now becomes 



while Eq. ( J4.8D becomes 



\ i j 
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These equations are both linear in y 2 . When y 2 is eliminated, the result is a factorizable cubic equation for x 

= x[x 2 (2 - 3pi) - x(4 + 2 Ps - 5pj) + 2(1 + Ps - Pl )] 



x(x — 1) 



x-2 



I+Ps-Pi 



2-2.pi 

Since pi = 1 — ps — Pm < 2/3, the physical solution is obviously 



(2-3pj). 



2 1 + ps - pj _ 2(2p 5 q 



2-3pj 



2-3p 7 



(4.34) 



(4.35) 



This leads to an expression for y 2 which is a quotient of somewhat complicated polynomials in pm, Pi, Ps- Those can 
be factorized, after some effort, leading to 



y 



(2 - <3 PI )(2 - PI )( PI - 2p s ) 
p 2 j{2p s +pi) 



which is positive if pi > 2p§. These results lead to the following expressions for «0j Po 

(2ps + Pm) (1 + Pm + Ps)(pi/2 - ps) 



pV> 



a 



p I {l-p I ) 2 +ps[2-{2-p I f] 

1 /2 

{pi - 2p s ){2ps + Pm +Ps )\ 

2-3 P i 

. pAps +Pm/2) 

Pi{l-Pi) 2 +Ps[2-{2-pi) 2 Y 



(4.36) 



(4.37) 



(4.38) 



Obviously, the cubic polynomial which appears in the denominators of oto, (3q is positive over the entire range < 
2ps < Pi where Eqs. ( [1.37] ) and (4.38) are applicable. Finally, Ao can again be found by first solving the following 
transcendental equation for the eccentricity e of the oblate spheroid 



1 + e' 



(e — arctane) =ps +Pm, 



and then using the relation between e and a, (3, A to get 



a 2 +(3 2 
a A 



- 1 



1/2 



aoAo 



1/2 



(4.39) 



(4.40) 
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These asymptotic results are valid for the range of constituent compositions defined by 



Ps < y < l- (4.41) 



When ps approaches pi/2 from below, then ao, fio, Ao all tend to 0, but at different rates: 

ao oc 



VPi ~ 2 VS, Ao S -22— ex y/ Pl - 2 PS , (4.42) 
1 + e z 



f3 <xpi~2p s . (4.43) 

It is worth noting that a, /3, and A are all independent of v. This indicates that the local electric field in the M 
constituent has considerable components that are perpendicular to B. It is also worth noting that, when pi — > 2/3, 
both ao and Ao diverge but ao/Ao — > 1, i.e., 

This indicates that, when the total volume fraction of conducting constituents j>m + Ps approaches its percolation 
threshold p c = 1/3, the same current flow paths in the M constituent are responsible for the leading contribution to 
the macroscopic response whatever the direction of the average current density (J) . 



2. Two -constituent M/I mixture 



An important special case within the nonsaturating regime is ps = 0, corresponding to a two-constituent M/I 
mixture. One can also get the results for such a mixture by setting ps = in Eqs. 
following equations: 



), (4.4), and (4.8) to get the 



X = u 



1 



(3 f Pi D 



M 



a 2 + f3 2 \pm Di 



PM ~ 

H 
1 + ff 2 



a f (3 > 2 

-A a -H ]=ia 



i + |)-^(" + l 



(4.45) 
(4.46) 

(4.47) 



The asymptotic behavior is obtained by making the ansatz of Eqs. (4.26). After a sequence of algebraic steps similar 
to those described above, the following results are obtained for the asymptotic linear slopes: 



00 



1 +PM 

2p 



on = 



Pi 



A I 



1 + PM 



2pm \ 3pM - 1 



1/2 



A = a 



(2p M jpi 
l + el 



where eo is given implicitly in terms of pm by the transcendental equation 



1 



Pm 



■ (eo — arctaneo), 



(4.48) 
(4.49) 

(4.50) 



which must be solved numerically. As in the M/I / S nonsaturating regime, ao, 0q, Ao are all independent of v. This 
indicates that the leading contribution to the macroscopic or bulk effective response is due to local currents that flow 
perpendicular to B in the M constituent. 

When pm approaches the SEMA percolation threshold value of 1/3 from above, then [3q — > 2 and eo — > 0. Therefore 
Xo/ao — ► 1 from above, and both ao and Ao diverge as 1/\/3pm — I. The fact that Ao/ao — > 1 again indicates that the 
same current flow paths in the M constituent are responsible for the leading contribution to the macroscopic response 
whatever the direction of the average current density (J) . 

Once again, in order to confirm the asympto tic b ehavior predicted analytically, we have solved the SEMA equations 
numerically. We assume V\ = u = 1 [see Eq. (3.1)], p% — oo, and consider a variety of values of p\ = pu above the 
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percolation threshold p c (equal to 1/3 in the SEMA in three dimensions). The resulting behavior of a and A is 
shown in Figs. ||(a) and|^(b); we also show the relative Hall resistivity (3/H in Fig. ||(c). Evidently, for any choice of 
Pm, ot and A rapidly approach their asymptotic linear dependence on |iJ|, as predicted by the asymptotic analysis. 
Furthermore, the slope increases as pm approaches the percolation threshold, again as predicted by the asymptotic 
results. The asymptotic linear dependence appears to be reached approximately when \H\ > 5, for all the values of 
Pm that were considered. For pm less than about 0.5, we had some difficulty, using our simple algorithm, in solving 
the SEMA equations numerically; by contrast, of course, the asymptotic analysis gives the slope, for any value of pm 
greater than p c — 1/3, without any difficulty. 



D. Scaling behavior near the transition point 



If w e compare the critical beha viors exhib ited by a, f3, A when \H\ ^> 1 and pi — > 2ps from above [Eqs. ( 1.42 ), 
( p3|) 1 and from below [Eqs. (fDJ) and flOcj)], we are led to anticipate a scaling behavior. That is, we should be able 



to describe the critical behavior in both regimes by using a scaling variable which is some power of H 2 {pi — 2ps) 3 /v 2 , 
and writing a, [3, A in terms of three (scaling) functions of that variable. We have found that the most convenient 
scaling variable for this purpose is 



\H\ 



v 



2/3 



Z=[!-L) ( PI - 2 PS ). (4.51) 



In the "critical region", i.e., when \H\ 3> 1 and \pj — 2ps\ -C 1, the bulk effective resistivity parameters can now be 
expressed as 



v 



F a (Z), (4.52) 

Pi - 2p s 

\= V — F X {Z), (4.53) 

pi - 2ps 



As usual, there are three important extreme regimes within the critical region, namely: Z < 0, \Z\ 3> 1 (Regime I, 
where 2ps > pi); Z > 0, \Z\ ^> 1 (Regime II, where 2ps < pi); and \Z\ -C 1 (Regime III, where we can have either 
2ps > pi, or 2ps < pi, or 2ps = pi). The behavior of the scaling functions F a (Z), F\(Z), Fp{Z) in Regimes I and II 
is dictated by the critical behaviors found earlier. Their behavior in Regime III is dictated by the requirement that 
this behavior must exactly cancel the divergences that would otherwise occur due to the vanishing 2ps — pi factor in 
the denominators of the above expressions. These considerations lead to the following forms for the scaling functions 
in the three regimes: 

-A Regime I 

F a {Z)={ A'Z 3 / 2 Regime II (4.55) 
A"Z Regime III 

—L Regime I 

F X {Z)={ L'Z 3 / 2 Regime II (4.56) 
L"Z Regime III 

B Regime I 

F () {Z)={ B'Z 3 Regime II (4.57) 
B"Z 2 Regime III 

where the primed and double primed and unprimed versions of A, B, L are positive dimcnsionlcss constants of order 1. 
Their values can be found, if necessary, by comparing the resulting expressions for a, [3, A with the detailed solutions 
of the SEMA equations. Qualitative plots of these scaling functions are shown in Fig. ||. 

As indicated above, these equations make a nontrivial prediction about the behavior of a, (3, A in Regime III, which 
was not worked out in the previous sections, namely: 

a = a \H\ 2 ' 3 , \ = X \H\ 2 / 3 , 13 S (3 \H \^ 3 sign(iJ). (4.58) 
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We now use these scaling forms as an ansatz for another asymptotic solution of the SEMA equations, assuming 
pi = 2ps, in order to check f or c onsistency of the scaling scheme developed here. 
Since A is very large, Eq. (4.3) entails 



n z =ps +PM = 1 - Pi = 1 - 2p s > ^, n x — Ps < T^j 



(4.59) 



where we have used the fact that pi = 2ps- From Eq. (4.3) we can also calculate the small correction to n z = 1 — 2ps 
in terms of Ao 



1 - 2p s - n z ^ 



2 Ps (l-3p s ) 



A |ff| 2 /3 l-2p s 



(4.60) 



In addition, we can write expressions for Dm, Di, and for their ratio Dm/Di, that go beyond leading order in 1/H, 
namely 



D M = (1 - n x f 



Di = (l- n x ) 2 + 



al\H\V3 

2«2 



a§|ffp/3' 



a/ 



Di 



^ 1 



1 - n ;7 



From Eq. (4.4) we get a relation between /3 and ao 



2 l-3p s 

Pa = a - ■ 

1 - ps 



(4.61) 
(4.62) 
(4.63) 

(4.64) 



Using this result, we find that Eq. (4.£) becomes an identity to leading order in 1/H. We therefore need to consider 
that equation in the next-to-leading order. The result is 



1 - 2ps - n z = 



a 2 Q 4p|(l-3p s ) 



\ H \2/Z (l_p 5 )3 ■ 

Comparison of this equation with Eq. ( 4.60] ) yields the following relation between ao and Ao 

QqAq = (1 -Ps) 3 
v 2p s {l-2ps)' 



(4.65) 



(4.66) 



Another relation between those two unknowns is obtained by first solving the transcendental equation for the eccen- 
tricity e of the oblate spheroid 



1 + e z 



(e — arctan e) = p s + pm , 



and then using the expression for e in terms of a, (3, A to get 

l + e» = ^ S 
qA 



Ao' 



(4.67) 



(4.68) 



Eqs. (4.64), (4.66), and (4.68) then provide a complete and consistent solution for ao, Po, Ao, obtained by jointly 
considering the scaling and SEMA equations. 

Again, it is worth noting that both a and A are proportional to is 1 / 3 , while (3 oc v 2 / 3 . This seems to be consistent 
with the fact that a and A are also proportional to |iJ| 2 / 3 , while (3 oc \H | 4 / 3 /iJ. These behaviors indicate that the 
local current flows both parallel and perpendicular to B, and that the three principal Ohmic conductivities of the M 
constituent all contribute to the macroscopic response when pi — 2ps- 
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V. SUMMARY AND DISCUSSION 



A striking result of the present work is that, in an M/I composite, both p^ and p^ are proportional to \H\ 
in the strong- field limit. Such linear magnetoresistance has long been a mysterious observed feature of transport in 
polycrystalline samples of even so-called "simple metals" .Ej This behavior has sometimes been ascribed to macroscopic 
inhomogeneities. But until now, only in the low-c oncen t ra tion limit has proof been given that such inhomogeneities 
could actually produce a linear magnetoresistance.liirEillia Here, we have shown that both p^ and p^ remain linear 
in \H\ even at higher concentrations of inclusions, provided those inclusions have strictly zero conductivity. This result 
may be relevant to a range of experimental systems. 

It is quite easy to understand, qualitatively, the behavior of an M/ S disordered mixture: Whatever the direction 
of the average current density (J), the local current density J(r) in the M constituent will flow mostly along B || z 
when H is very large, and this tendency will become more and more pronounced with increasing \H\. Components 
of J(r) that are perpendicular to B will have finite values only inside the S constituent, while in the M constituent 
they will tend to as l/H 2 . This results in a current flow pattern that saturates when \H\ ^> 1; therefore, all the 
Ohmic resistivities will also saturate. The Hall resistivity will be very small — of order l/H — because the local Hall 
field ~Eh will also be of that order. 

The behavio r of M\/Mi mixtures is more difficult to understand qualitatively. Our interpretation of the (successful) 



ansatz of Eq. (3J3) is that when (J) || B || z, the current distribution saturates for large values of \H\. But when 
(J) _L B, then that distribution continues to evolve with increasing |iJ|, with the current distortions increasing 
asymptotically as If/] 1 / 3 . 

The macroscopic response of M/I mixtures is qualitatively similar to the behavior of an isolated insulating inclusion 
embedded in an M host, as shown many years ago by an exact solution for such an inclusion of spherical or cylindrical 
shapeH The fact that this kind of behavior persists down to the percolation threshold indicates that, despite the 
interactions between current distortions produced by different inclusions, the local current distribution never saturates 
as \H | increases. The results we got would require that the current distortions increase as iTf] 1 / 2 for large H. However, 
the fact that the coefficients ao and Ao diverge as pm — ► Pc probably signals that those distortions increase even more 
rapidly than If/] 1 / 2 precisely at the percolation threshold. 

The present results can also be compared with some simulations performed previously an a discrete network model, 
where finite size Lx Lx L samples were considered precisely at the percolation threshold.E 2 ] Those simulations showed 
that, when L is much less than a "magnetic correlation length" oc |i/| ' 46 , the Ohmic resistivities saturate at 
a value proportional to L . However, in the opposite limit L 3> (,h , the Ohmic resistivities continue to increase 
as |ff| 21 i 2 - 2 . This result is consistent, within the error bars, with 7J 2 L*/ y , where t = 2.0 and v = 0.89 are the 
usual percolation critical exponents for scalar Ohmic conductivity (<7 e cx Ap*, Ap = pm — p c ) and for the percolation 
correlation length (£ p oc Ap~ u ) — see Ref. 0. Using finite size scaling to replace the system size L by £ p in these 
results, we are led to expect that 

/ Ap~ 4 for Ap » |iJ|-°- 52 or £ p « fr, 
a > Aoc iff2 for Ap« |iJ|-°- 52 or^»^. [ > 



Obviously, this behavior differs in a number of ways from what was found in Section IV C 2 using SEMA. It is 
not surprising that the critical exponents predicted by SEMA are quantitatively incorrect. However, the fact that, 
according to the simulations, the magnetoresistivity saturates as \H\ — > oo, is qualitatively at odds with the predictions 
of SEMA. Clearly, this qualitative discrepancy needs to be examined further. We conjecture that it has to do with 
the existence of three diverging lengths in this problem, namely £ p , £#, and L. 

The qualitative situations described above continue to be applicable also in the case of the three-constituent M/I/S 
mixtures. The presence or absence of system-spanning (i.e., infinite) current flow paths, which are parallel to B inside 
the M constituent, now depends on the relative amounts of I and S inclusions. The critical points or transition points 
pi = 2ps can be expected to occur when such paths first appear with increasing ps- 

Obviously, this transition is a kind of percolation process. In fact, we believe it is a physical realization of "anisotropic 
percolation" . This kind of percolation— svas first considered many years ago in the context of a two-dimensional, 
randomly diluted, square bond network. Ej The geometrical propcrties-cif-anisotropic percolation in hypercubic random 
bond networks of arbitrary dimension were also studied extensively.L§E3 

As originally defined, the anisotropic percolation problem is characterized by different occupation probabilities for 
bonds aligned along different principal axes of the network. This situation is not easily implemented in a continuum 
percolating system, because that would require an anisotropic constituent where the principal conductivities have 
ratios that are extremely different from 1. Also, the principal axes of different grains of that constituent would have 
to be identically oriented. While this may be difficult to achieve in conducting materials when B = 0, such extreme 
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ratios and identical orientations can easily be attained, even in an isotropic conductor, just by applying a magnetic 
field such that \H\ > 1 or v< 1. In that case, if we identify volume fractions with the bond occupation probabilities, 
then we can say that, for \H\ ^> 1 or v <C 1, the regions that are highly conducting along B occupy a fraction pu + ps 
of the total volume. By contrast, the regions that are highly conducting perpendicular to B occupy the smaller 
volume fraction p$. Our SEMA result, for the transition points between saturating and non-saturating regimes of 
magnetoresistance, can then be interpreted as follows 

= 1 - p x - P y - Pz = 1 - Ps ~ PS - (pm + Ps) 
= Pi- 2ps, (5.2) 

where p z = pm + Ps represents the bond occupation probability along B while p x = p y = ps represent the bond 
occupation probabilities in the directions perpendicular to B. 

The identification of the critical points in the M/I/S magnetoresistive response with an anisotropic percolation 
threshold needs to be verified by a more accurate treatment of the bulk effective magnetoresistance. In particular, it 
needs to be determined whether the ever-present Hall conductivity is an irrelevant perturbation, or whether it in fact 
destabilizes the usual percolation fixed point, or alters the critical behavior associated with it. In any case, we expect 
that such a treatment will yield different values for the critical exponents and different forms for the scaling functions 
and scaling variable than those obtained here using SEMA. 

Experimental studies of the line of magnetoresistive critical points can be conducted on samples made by using a 
semiconductor host with a large Hall mobility fi as the M constituent, in which two types of inclusions are randomly 
embedded: highly insulating inclusions (e.g., voids) as the / constituent, and either highly conducting normal metallic 
inclusions or superconducting inclusions as the S constituent. 
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FIG. 1. Results of numerical solution of the SEMA Eqs. ( 3.3 )— ( 3.5) for an M1/M2 mixture, (a) Plot of the relative bulk 
effective transverse magnetoresistivity a — 1 vs. |i?| 2//3 for a three-dimensional composite of two free-electron-metal constituents, 
present in equal amounts pi — p2 = 1/2, which have the same ohmic resistivity, but Hall resistivities that differ by a factor 
two [i.e., pi = p2 and Hi = H , H2 = 2H, in the notation of Eq. (3.1)]. (b) Same as (a), except that we plot the relative bulk 
effective Hall resistivity (3(H) /H vs. \H\ 2/3 [cf. Eq. (E|)]. 
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FIG. 2. Results of numerical solution of the SEMA Eqs. ( 4.45 )-( 4.47) for an M/I mixture, (a) Plot of the relative bulk 
effective transverse magnetoresistivity a — 1 vs. |_ff| for a three-dimensional two-constituent composite of free-electron metal and 
perfect insulator for various values of the metal volume fraction p = pu- (b) Same as (a), except that we plot the relative bulk 
effective longitudinal magnetoresistivity A — 1. (c) Same as (a), except that we plot the relative bulk effective Hall resistivity 
/3(H)/ H vs. \H\ only at high fields (\H\ > 5). 




FIG. 3. Qualitative plots of the scaling functions obtained using SEMA. The left plot shows F a (Z) (solid line) and F\(Z) 
(dashed line), which are similar but not identical, with \F a (Z)\ > |Fa(Z)|, while the right plot shows Fp(Z), which is never 
negative. 
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